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The calculation of partial decay rates in B ^ Xuf-V decays at next-to-next-to-leading order (NNLO) in as and to 
leading order in l/mj, is described. New results for the hard function are combined with known results for the jet 
function and shape-function moments in a numerical analysis which explores the impact of the NNLO corrections on 
partial decay rates and the determination of IVutl- 



1. INTRODUCTION 



The inclusive decay B Xuiv is of much interest because of its potential to constrain the CKM element \Vub\- 
Due to experimental cuts required to suppress charm background, measurements of this decay are available only 
in the shape-function region, where the hadronic final state is collimated into a single jet carrying a large energy 
on the order of and a moderate invariant mass squared on the order of mftAqcD- The theory challenge is to 
^ calculate partial decay rates in the presence of these cuts, where a local operator product expansion is insufficient. A 
• systematic treatment relies on a non-local operator product expansion whose end result can be formulated in terms 
of a factorization theorem. Different approaches to this factorization have been put forth in the literature, going 
ji^, under the names of BLNP [l|, Q , GGOU 3| , and the dressed-gluon exponentiation ^ . The purpose of this talk is to 
^ ' describe the elements that go into the BLNP formalism at next-to-next-to-leading order (NNLO) in and leading 
order in 1/m;,. 

The main point of discussion is the factorization formula for an arbitrary decay distribution restricted to the 
kinematics of the shape-function region 

^; dr^H-j<i^s + o(^^^), (1) 

C^l ' where the symbol (8" denotes a convolution. The factorization formula contains a hard function H, which is related 



to physics at the hard scale rrib, a jet function J, which is related to physics at the intermediate scale mfeAqcD, and a 
non-perturbative shape function S, describing the internal soft dynamics of the B meson 0, 0]- This factorized form 
was originally derived in 0, B| using diagrammatic techniques, and was rederived in the framework of soft-coUinear 
j-jl ' effective theory (SCET) in [if, @|. The hard and jet functions to NLO in perturbation theory have been known for 
■ some time Q, and the NNLO jet function was obtained in [lO|. Very recently, the NNLO contributions to the 
, hard function have also been calculated [ll|, [l2, 13 1, thus completing the perturbative corrections to the factorization 



Ci formula to this order. In Section[2]we review the calculation of H and J to NNLO, along with the form and solutions 
of their renormalization-group (RG) evolution equations. Then, in Section [3l we give preliminary results illustrating 
the numerical impact of the NNLO corrections on the P+ spectrum, comment on the relevance for the determination 
of \Vub\, and make some concluding remarks. 



2. THE HARD AND JET FUNCTIONS AT NNLO 

In this section we describe the calculation of the hard and jet functions at NNLO, as well as the form and solution 
of their RG evolution equations. 

The hard function arises when integrating out fluctuations at the scale nib by matching QCD onto SCET. This 
matching can be done at the level of the b ^ u transition current, and to leading order in l/rrib takes the form 

3 

e-™''"-"M(a:)7^(l - ro)b{x) = J2 dsC^{s)x{x + sn)rfH(x-), (2) 

i—1 
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where we have followed the SCET conventions of 1], and the are a set of three Dirac structures. In practice, the 
matching calculation is carried out in momentum space and yields results for the Fourier-transformed coefficients, 
which read 



Ci{n-p) 



(3) 



The hard function H is derived from the matrix of coefficients Hij — CiCj. 

To obtain the matching coefficients Ci requires to calculate UV-renormalized matrix elements in full QCD and 
SCET. The calculation is simplest when the external states are chosen as on-shell quarks and both UV and IR 
divergences are regulated in dimensional regularization in d = 4 — 2e dimensions. In that case the loop corrections 
to the SCET matrix elements are given by scaleless integrals and vanish, so that the result is just its tree-level value 
multiplied by renormalization factors from operator and wave-function renormalization. The main challenge is to 
calculate the QCD result, which is written in terms of three Dirac structures multiplied by scalar form factors Di. 
To obtain these scalar form factors at NNLO requires to calculate the two-loop corrections to the b u current in 
QCD. This task has been completed in [ll|, [l2, 13|. One then turns the results for the scalar form factors Di into 
those for the SCET Wilson coefficients by evaluating the matching conditions 



Ci{mb,n ■ p, fi) = liia Zj^{e, nib, n ■ p, fi)Di{e,mb,n ■ p, fi) 



(4) 



where Zj is a current renormalization factor. Details of the matching calculation and results for the Wilson coefficients 
have been presented in l3, 

In the BLNP approach, the hard function is evaluated at a scale Hi ~ 1.5 GeV and logarithms of the ratio fJ.h/^J■^ are 
treated as large. One can resum these logarithms by deriving and solving the RG equation for the hard coefficients 
Ci, or equivalently the matrix of coefficients Hij. The evolution equation has the form (see, e.g. [l|) 



dln/z 



Hij (n • p, /i) = 2 



7' (as) -f rcusp(as)ln 



n ■ p 



H,j{n-p,fi) , 



and its solution is 



exp 



25(^;i,/i) - 2ar(Ai/i,Ai) In 



2a-y/(/i/i,/i) 



Hij{n-p, y.h) , 



(5) 



(6) 



where we have defined y = n ■ p/nib. Explicit results for the Sudakov factor S and the anomalous exponents ar, a^' 
can be read off from A consistent treatment at NNLO in as requires the matching coefficients at two loops, the 
anomalous dimension 7' at three loops, and the cusp anomalous dimension at four loops. The anomalous dimensions 
are both known to one loop lower, which adds a small uncertainty to the analysis. 

We now turn to the jet function J, which arises when integrating out fluctuations at the intermediate scale mbAqcD 
by matching SCET onto HQET. The calculation of the jet function can be recast into the evaluation of the imaginary 
part of a certain vacuum matrix element in QCD, and was calculated to NNLO in [l3|. In light-cone gauge, this 
matrix element is the quark propagator. The jet function depends on a renormalization scale /i, which is usually 
kept fixed at /ii = 1.5 GeV in the BLNP approach. However, it is possible to derive and solve the RG equation for 
J, which allows one to also vary the renormalization scale at which the jet function is evaluated. The starting point 
is the exact integro-differential evolution equation derived in [10| : 



cgJ(p^M) 
dln/i 



2F, 



cusp In ^ -f 27'' 



J{p\^JL)-2T. 



cusp 



dp 



p- ~p 



/2 



(7) 



A novel solution to the evolution equation ([7]) was presented in 15|. The method is based on the observation that 
the Laplace transformed jet function 



J 



dp^ e 



1 



(8) 
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Figure 1: Partial rates with a cut P+ < M^/Ms as a function of fin at LO, NLO, and NNLO. 



obeys a local evolution equation analogous to ([5]). Solving as before and inverting the Laplace transform gives the 
resummed J in momentum space. Alternatively, it is always possible to write partial decay rates in terms of an 
integral over the jet function, defined as 

j(^n^^^= dp^J{j>^p:}. (9) 
Using the solution of the RG equation, one can write a compact expression for the resummed integrated jet function 

where expressions for the RG factors a^j and rj can be found in [l^. Equations ^ and pO|l allow one to evaluate 
integrals of the product iJ • J at an arbitrary scale /x/ at which the shape function is renormalized, and study the 
dependence of partial decay rates under independent variations of the matching scales fih and fii, under which they 
are formally independent. Detailed numerical results shall be given in 14 1. 



3. NUMERICAL RESULTS AND CONCLUSIONS 

In this section we briefly explore the numerical impact of the NNLO corrections, taking as an example the partial 
rate with a cut P+ < A M^/Mb- This partial rate can be written as 

r„(P+ < A) = Qg^3 dujS{uj,Hf) dq5q dyy ^Hu{y,Hf)j l^n , ^/ j (11) 

where A^ — min(A, Mg — q), and the hard and integrated jet functions are to be evaluated in their resummed forms 
© and (fTU]) . The function is related to certain combinations of the Hij and can be deduced from . To evaluate 
(|lip requires a model for the shape function. In what follows, we use the two-parameter exponential model 

S{Lb) = M{b, A)(2-''-iexp (^-^^ ■ (12) 

It is possible to put model-independent constraints on the normalization factor J\f and the parameters b and A by 
studying shape- function moments [l[. These moments are defined as 

M]y{ujo, fJ.i) = / dujuj^ S{ui,Hi) , (13) 
Jo 

and can be obtained as an expansion in local HQET operators as long as wq 3> Aqcd ■ The zeroth moment fixes the 
normalization of the shape function and was determined to NNLO in [l6l|, and can be used along with the first and 



3 



34* International Conference on High Energy Physics, Philadelphia, 2008 

second moment to define the parameters mi, — Mb — A and /i^ in the shape- function scheme. Numerical values for the 
HQET parameters nii, and /xj in the shape-function scheme at NNLO can be determined from those obtained from 
global fits in other schemes using the two-loop conversion relations derived in [l7| . One then tunes the parameters 
b and A such that the moments of the model function (jl2[) give appropriate results for mf, and /i^. In our analysis 
here, we tune the parameters such that mt, = 4.61 GeV and /ij — 0.2 GeV^ at NNLO in the shape-function scheme. 

In Figure [TJ we show results for the partial rate with a cut P+ < A = A/^/A/s, for the choice fii — fif — 1.5 GeV, 
as a function of the hard matching scale iih- From the figure, one sees that the dependence on the renormalization 
scale fifi is reduced when going to higher orders in perturbation theory, and also that there is a fairly large downward 
shift in the central value between NLO and NNLO. This would tend to raise the value of \Vub\ compared to the result 
deduced from the NLO calculation in the BLNP approach. However, we must stress that these numerical results are 
preliminary, and will be finalized in [14| . along with uncertainties associated with variations of the matching scale 
Hi. Results for other partial rates, such as those with cuts on the hadronic invariant mass and lepton energy, will 
also be presented. 

It is worth emphasizing that the partial decay rates are rather sensitive to the numerical value of mf,. If, as 
suggested in [l8|, the B Xs^ moments are excluded from global fits on the grounds that measurements of the 
photon energy spectrum are made in kinematic regions where shape- function effects are important, then mt, tends to 
come out closer to 4.71 GeV. Changing rrib by such an amount raises the partial decay rates upwards by approximately 
20%. A precise determination of \Vub\ from inclusive decays will thus require that this point be settled. 
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